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Following the development of a scheme to bosonize and debosonize consistently [N. Shah and C.J. 

Bolech, Phys. Rev B 93, 085440 (2016); arXiv:1508.03078], we present in detail the Toulouse-point 
analytic solution of the two-lead Kondo junction model. The existence and location of the solvable 
point is not modified, but the calculational methodology and the final expressions for observable 
quantities change markedly as compared to the existent results. This solvable point is one of the 
remarkably few exact results for nonequilibrium transport in correlated systems. It yields relatively 
simple analytical expressions for the current in the full range of temperature, magnetic field, and 
voltage. It also shows precisely, within the limitations of the Toulouse fine-tuning, how the transport 
evolves depending on the relative strengths of interlead and intralead Kondo exchange couplings 
ranging from weak to strong. Thus its improved understanding is an important stepping stone for 
future research. 


I. INTRODUCTION 

In the first part of this series [1] , we introduced a con¬ 
sistent prescription in order to be able to bosonize, make 
transformations, and debosonize consistently in the pres¬ 
ence of “active local impurities or boundaries”, which we 
called the consistent bosonization-debosonization (BdB) 
program. In this paper we explore the implications of 
this formalism for the important case of quantum impu¬ 
rity problems. 

Just over fifty years ago, in 1964, Kondo showed that a 
mysterious finite-temperature minimum in the resistivity 
of metals was due to the contributions from dilute mag¬ 
netic impurities present in the samples [2]. This marked 
the start of the study of the so-called Kondo problem, 
which is one of the pillars of modern condensed-matter 
theory [3]. On the technical side, the reason for this is 
that the problem of a single magnetic impurity in a metal 
is one of the very first examples of an asymptotically free 
theory [4]. When the system is below a certain energy 
scale known as the Kondo temperature, standard per¬ 
turbation theory fails and one needs to resort to more 
sophisticated theoretical tools (many of which were ac¬ 
tually first developed studying this problem [5]). These 
range from the exact to the versatile, or from the Bethe 
ansatz [6-10] to auxiliary-particle perturbation methods 
[11-15]. On the application side, the relevance of the 
Kondo problem extends nowadays well beyond the origi¬ 
nal system of impurities in metals. The so-called Kondo 
lattice is the central model in the study of heavy fermions 
[4] and, even more generally, a formalism know as dynam¬ 
ical mean field theory is based on the mapping of any 
complicated tight-binding model of a material to an ef¬ 
fective quantum impurity problem [16]. Moreover, since 
the last two decades, as the study of artificial mesoscopic 
systems reached the nanoscale, the Kondo problem can 
show up in all sorts of electronic devices, most notably in 
semiconductor quantum dots [17, 18] but also in molec¬ 
ular electronics [19], etc. 

The experiments with artificial nanostructures brought 


in an additional layer of complication to the theory of 
the Kondo effect. Most of the typical experiments in¬ 
volve transport measurements in nonequilibrium condi¬ 
tions, while our best theoretical tools to describe systems 
out of equilibrium are perturbative. Since the Kondo ef¬ 
fect is nonperturbative, a lot of theoretical activity con¬ 
tinues to ensue, to the point that it is fair to say that a 
deep understanding of correlated systems out of equilib¬ 
rium is still work in progress. 

Early on during the nanoscale revolution in mesoscop- 
ics, the pioneering theoretical work of Schiller and Her- 
shfield (S&H) provided the first (and still now one of the 
few) exact solution of a nonequilibrium strongly corre¬ 
lated quantum problem [20]. A few years before them, 
Emery and Kivelson (E&K) found a solvable point (called 
a Toulouse point) for the two-channel Kondo model [21] 
(in which, besides spin, the band electrons have one more 
two-valued discrete degree of freedom [22]). S&H realized 
that they could adapt it to the case of a Kondo impurity 
interacting with two separate leads. They thus found a 
mapping that allows for the calculation of transport in 
a problem that has strong correlations due to the ex¬ 
change interaction between the electrons and the impu¬ 
rity spin. We will revisit the central aspects of their work 
below. We shall find that while the key insights of S&H 
about the existence of a nonequilibrium Toulouse point 
remain valid, our generic consistent-debosonization pro¬ 
cedure shows that the actual observables being calculated 
are substantially modified and yield more physically con¬ 
sistent results. 


II. MODEL OF A KONDO JUNCTION 

We shall be interested in a class of systems in which the 
transport between two leads or terminals happens across 
a microscopic region in which the fermionic degrees of 
freedom are such that the region has a total magnetic 
moment that remains unscreened (and we will focus on 
the case of spin-1/2). This situation is typical of nano- 
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FIG. 1. Schematic depiction of the setting in which two Fermi 
seas kept at different chemical potentials with their difference 
given by eV = pn — fJ-B. are connected via quantum tunneling 
across the potential barrier that separates them. The situa¬ 
tion when the barrier region does not allow for internal states 
was discussed before [1]. Here there exists a many-body state 
trapped by the barrier with an unscreened total spin 1/2 that 
we refer to as the impurity (depicted by an arrowhead). Due 
to a strong Coulomb blockade, the impurity interacts with 
the electrons in the leads via exchange processes only. For 
the purpose of the figure we generically denote by Jra the 
intra-lead exchange terms (either parallel or perpendicular) 
and by Jer the inter-lead cotunneling exchange processes that 
can give rise to a current (a similar notation will be introduced 
in the text later). 

scale quantum dots with strong Coulomb blockade and 
has been an experimental reality since the late 1990s [17, 
18]. At low temperatures, remarkably, the system enters 
the so-called Kondo regime and is able to conduct despite 
the Coulomb blockade. 

We are thus interested in the low-temperature char¬ 
acteristics of nonlinear transport across a quantum dot 
in the Kondo regime. We will model the system with 
a two-lead version of the Kondo model (see Fig. 1). In 
analogy to the equilibrium case, this model can be de¬ 
rived from a more microscopic Anderson-type model via a 
(time-dependent) Schrieffer-Wolff transformation [5, 23]. 
In Hamiltonian language, the model is given by 

H = 'y ] ( f THPi dx + A Afield , (la) 


(7,i '' ' 

t) (-ivpdx) [x, t) , (Ib) 

= Ju' (fV'L (0, t) (0, t)) , (Ic) 

= Jw SLb (0, t) (0, <)) , (Id) 

Tffield — ^‘^imp ; (1®) 


where {x, t) are chiral fermions in the Heisenberg rep¬ 
resentation that are obtained after unfolding the two 
leads in the usual way [1, 24]. We adopt the notation 
cr = {i,t} = {-1)+1} and I = {L,R} = {-1,-|-I} for 
the spin and lead index, respectively. The bar notation 
denotes a sign change—e.g., tf = —a. The impurity is 
described by and S'F^p = S^^p + iaSf^p. We assume 
now that at a much earlier time the connection between 
the two leads was established and that there is a battery 


keeping a constant chemical-potential difference between 
the two leads (what we call a Landauer-type configura¬ 
tion [25, 26]). Under these conditions the system will be 
in a nonequilibrium steady state [27]. Let us call /if the 
chemical potential of lead £, such that /ip — /^r = eU, 
with V the voltage drop across the junction. The in¬ 
formation about these chemical potentials will enter into 
the distribution functions for each lead. 


III. BOSONIZATION-DEBOSONIZATION 
APPROACH 

To set the stage for the bosonization of the model, we 
start by gauging away the chemical-potential difference. 
For that, let us first switch to Lagrangian language in 
which the system is described by 

= V'L (a;, t) (idt) (x, t) - 

= V’L {x, t) (idt + ivpd^) {x, t) , (2a) 

Lk = -Hk = -Jh' ^Lp (0, t) i/.f, (0, t)) , (2b) 

= -Hk = -Ju' STnip QV'L (0,i) (0,t)^ . (2c) 

We can now make the following (gauge) field transfor¬ 
mation ipa-e (x, t) = (x, t). The important point 

is that now the distribution functions do not contain 
information about the chemical potentials any longer 
(cf. with the discussion for the case of a simple junc¬ 
tion [1]). Next we subtract the vev (vacuum expectation 
value), which for a noninteracting problem is equivalent 
to factoring out the fast oscillations in each lead accord¬ 
ing to ip^e{x,t) = (x,t), with fcf = /if/up for 

this linear-dispersion case. So we are naturally, thanks 
to the linear dispersion, lead to the normal-ordered for¬ 
mulation of the problem: 

(^> ^) {x, t) : , (3a) 

Lk = -Jh SLb (I: (0. t) (0, t) : ) - 

- 5f,„p (0, t) (0, t)) , (3b) 

= ~Ju '5'imp ^2^^^ ~ 

- 0 (0, ■ (3c) 

Notice that in the case of the parallel intralead im¬ 
purity terms we wrote them also as normal ordered 
(since the vev’s of the two spin projections cancel each 
other due to the ct/ 2 factor), which is customary in 
other approaches such as boundary conformal field the¬ 
ory (BCFT), whereas the interlead impurity terms nat¬ 
urally remain non-normal-ordered. At this point we lost 
the information about the absolute energy reference but 
we still have the information about the potential drop 
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encoded in the time-dependent phase of the tunneling 
term. 


A. Bosonization and Initial Mappings 

The hrst part of the BdB program starts by bosoniz- 
ing, of course. In order to do that, we go back to the 
Hamiltonian formulation of the problem, 


consistent BdB treatment of the simple junction problem 
[1]; notice as well that we are also introducing the 1/2 
factors associated with consistent boundary conditions 
(CBCs) [1]. We get 


m = 


J-i 
J B. 


2'Ka 2 


ncu-j 

'^imp 




AieVt 




2'Ka 


SZ 


( 




imp 


) , ( 8 ) 


(2^: i) {-iwdx) Ae {x, t) : , (4a) 

Hk = Ju (^2 ■ 

+ , (4b) 

= Ju -S'imp (0> t) (0^ + 

+ '8'imp (0> 0 ; (4c) 

and we bosonize according to H/, with taken as 

the interaction terms. We follow the same standard 
bosonization prescription as we did for the junction prob¬ 
lem [1], '^cre{x,t) = —J=Fcre (t) and in terms 

of the bosons the Hamiltonian density for the leads takes 
the usual form (the Klein factors, F^^, drop out from 
these terms) 

n° = = : [d^(j)ae{x,t)f : . (5) 

£ (T=t,i 

^=L,R 

Using the same standard physically-motivated rotated 
boson basis, (fig-e = i4>c + <x4>s + ^'Pi + o’^/'si) /2, as we did 
for the simple junction [1], the non-interacting Hamilto¬ 
nian density retains its quadratic form 

= ^ . ( 6 ) 


Let us postpone the discussion of and proceed to 
bosonize the “perpendicular” part of the Kondo term. 
The first part is the intralead one, or lead-nonmixing, 
and it is present also in the standard two-channel Kondo 
model, while the second part is interlead, or lead-mixing, 
and is responsible for transport as can already be seen 
from the voltage dependence (we keep the time and space 
dependence of the bosonic fields implicit for the sake of 
brevity): 


m = 


j, 


u 


27ra 


SL 


imp 


_pt i4>s 

2 


'F^ie 


icficj 


J-k 

^ileVt qcr 

27:a “ 


i pt p -p-i4>a 


(7) 


We now change to the rotated boson basis and being 
careful of not combining vertex operators with opposite 
signs we introduce h factors in the same manner as in our 


where almost all possible fi factors, i.e.. 



(9a) 


(9b) 


(9c) 


(9d) 


appear except for the ones from the spin sector. 

The same as in the case of the simple barrier junc¬ 
tion, we do not expect the Klein factors to modify the 
physics. We treat them as we did in that case [28, 29] by 
identifying relations between different bilinears of origi¬ 
nal and new Klein factors and fixing the four arbitrary 
phases; see Eqs. (16a)-(16d) from Ref. 1. The rest of the 
Klein-factor relations can be derived from these [30]. In 
particular, for the intralead terms, we need 

= fIf! , (10a) 

= FsiF^ , (10b) 

fIrF^k = FsFsi , (10c) 

FlFrL=FsFl,, (lOd) 

whereas for the interlead terms we need 

fI^F^^ = F}fI, (lOe) 

F^F^r = FiFj , (lOf) 

T^^l^TR = FiFs , (lOg) 

fIrF^l = E.F/ . (lOh) 


Replacing the Klein factors and expanding explicitly 
one arrives at 


J- _ Jrr Anf 


m = 


2'Ka 2 


Jfi Rch) 


27ra 2 


[SYFsFsie-^^^e-'‘^^^ + 


■ 9+ p -i 




-ieVt^KL o+ 


2Ka 2 




r^uFi^pF.FA -“e" 
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ieVt 

2'Ka 2 
+ fl 




i^iFte*^=e-**5rp) . (11) 


It is easy to notice that the impurity spin and the lead 
degrees of freedom associated with the spin sector always 
appear together in the combination S^pFse~^'^‘ and its 
Hermitian conjugate. 


B. Toulouse Limit and Completion of the Square 


It is natural to describe the strong-coupling limit be¬ 
tween the impurity and the electrons by looking for a 
transformation that binds those two degrees of freedom 
together as a single one. As shown by E&K, that is 
achieved by the following transformation (the boson field 
is evaluated at the position of the impurity): 

U = . (12) 


It follows by simple algebra that U is unitary and com¬ 
mutes with all the Klein factors and vertex operators 
(notice that one defines U so that no point splitting 
is involved when applying it). Using the spin alge¬ 


bra, 


C± CZ 

^imp ’ ^imp 


= tS- 


imp ’ 


and the Baker-CampbelU 

Hausdorff (BCH) formula, e~^Ae^ = A -\- [A, B] + 
^ [[A, i?], B] -I- ..., we find 


C/5± pl/t = , (13a) 

C/5f^pI/t = 5f^p . (13b) 


Therefore, the perpendicular Kondo term in the Hamil¬ 
tonian transforms as = UH^U\ which with the sim¬ 
plifying choice of 7 s = 1 to absorb the spin-sector vertex 
into the impurity, and further defining (so 

that d'^d = + 1/2), takes the form 


H/p = 




27ra 2 
LL ^cflT 


si + FlF'f>‘‘d 




(d/Fgie 


+ Fgie-^*^^d 


27ra 2 

'^R^L 
2'Ka 2 

- (nt, . (14) 

Note that no nf factors appeared in these terms; they 
would appear in the inter lead cotunneling terms of Fl^, 
but we follow S&H and set the coupling constants of those 
terms to zero as part of the definition of the Toulouse 
limit. So we need to discuss the intralead part of and 
the kinetic terms. Let us examine further the effects of 
the E&K transformation. The transformation of a boson 
derivative is given by 

Ud(j)s (x) = d(j)g (x) - ijsSiljp [d(j)s (x ), (/)g (0)] -k ... 

= d(/)g (x) - 27r7^S'f„p(5 (x) , (15) 


where we used Ae^ = e-^ (A + [A, H] -k ...), which fol¬ 
lows immediately from the BCH formula and the equal¬ 
time commutator (x ), (y)] = 2Ki6 (x — y) 

In terms of the corresponding fermions, via debosoniza- 
tion, this shift translates into a change of boundary con¬ 
ditions and thus U is sometimes called a “boundary- 
condition changing operator” [31, 32]. When 7s = 1 this 
gives (up to a conventional sign) an Abelian version of 
the shift that Affleck and Ludwig use to “complete the 
square” and absorb the impurity in a redefinition of the 
“spin density” at the infrared fixed point [24]. On the 
one hand, for the spin-sector part of the kinetic energy, 
we use this shift and obtain 

= y {d(t>s {x) - 27r7sS'f^p5 {x)f = (16) 

= ^ (*)]^ - (0) -k llv'p'^5 (0) . 

On the other hand, for the parallel Kondo terms we get 

= y [d</s (0) - 2K^gS/^p5 (0)] + 

+ ^ ('^RR ~ '^ll) , (17) 

where Jf^g = (J^r + J£l) /2. Eor the Toulouse limit 
one considers the symmetric case, J^p. = J£l, and sets 
—>■ 47ruF7s = 'Ikv'p- Combining these two results 
and disregarding constant energy shifts, one has 

U {Hl^g + H^) I/t = ^ [a^, (a,)]2 . (18) 

In summary, all the parallel Kondo terms were either set 
to zero or absorbed into the kinetic term and dropped 
out from the problem. 

Finally, the local-field term is not affected by the trans¬ 
formation procedure and is written as 

Ffleid = -h = -h {dU - 1/2) . (19) 

C. Debosonization 

The kinetic terms are easily written back in terms of 
fermions, becoming similar to the original kinetic terms 
of the model. The only nontrivial part of the Hamiltonian 
that we need to debosonize and discuss carefully is the 
“perpendicular” Kondo terms. The final result for the 
lead-symmetric case reads 

Hk = (d^'ipsi (0) -k V'L (0) d) + 

+ (0) + tpsl (0) d) - 

- (^d^i (0) + V'/ (0) d^ - 

- der ”y^ (W (0) dU d'i/l (0)) , (20) 
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where (cf. Fig. 1) 


as (notice = 0) 


Jra, = Ju /27ra and Jer = ; (21) 


and later we will compactly denote 


J± = 


Jra^c^l 
4:VF 


and 


^ 4wf 


( 22 ) 


The fi^’s defined in Eq. (9) are now viewed as square 
roots of local fermionic densities at the site of the im¬ 
purity. While writing Eq. (20), we have already gauged 
away the time dependence from the coupling constants 
so that = —{eV) and all other chemical potentials 
are zero. This is achieved by using the transformation 
ijji {x,t) = {x,t) where '0^ = 

are the debosonized fields, as was done for the simple 
junction [1]. 


IV. TRANSPORT CALCULATIONS 

In order to solve for the transport characteristics, we 
derive an expression for the current according to / = 
dt^ =i [H, = i [H^,N^=i\ , which gives / = (^/^ 


I = - 


. r IT'c 

% U or 

2 L 




(23) 




and the problem reduces to finding those matrix ele¬ 
ments. 


A. Conventional Approach 

In the conventional BdB program, the boson exponen¬ 
tials that result after changing basis are freely recom¬ 
bined and as a result they simply disappear. This is 
equivalent, in the expressions above which already incor¬ 
porate CBCs that are conventionally not discussed, to 
replacing nj —)■ I everywhere [cf. Eq. (9)]. 

To calculate the necessary Green’s function elements, 
using the same local action scheme and principal-value 
regularization (cf. Ref. 33) as we did for the case of 
the junction [I], we adopt the following Keldysh-Nambu 
spinor basis (with the frequencies restricted to the posi¬ 
tive semiaxis only in order to avoid double counting), 


4-(w) = (?/>; (w) {uj) V’l (w) ^/"^(w) d (w) d+{uj) (w) d^+(w) '(/’«; (w) (w) i/jh (w) ■ 


Let us define Si, (w) = tanh and Sjy (w) = tanh 2 t^" ? where Temp is the temperature that is taken to be uniform 

(and that we will set to zero for the most part). The local inverse Green’s function for the junction, G~^ (cu) /2vf, is 
thus given by the following matrix: 


isi 

-isi + i 

0 

0 

T+ 

0 

T_ 

0 

0 

0 

0 

0 

-isi - i 

isi 

0 

0 

0 

-T+ 

0 

-r_ 

0 

0 

0 

0 

0 

0 

isi 

—isi + i 

-T_ 

0 

-T+ 

0 

0 

0 

0 

0 

0 

0 

-isi - i 

isi 

0 

T_ 

0 

T+ 

0 

0 

0 

0 

T+ 

0 

-T_ 

0 

uj -\- h 

0 

0 

0 

-J+ 

0 

-J_ 

0 

0 

-T+ 

0 

r_ 

0 

—U! — h 

0 

0 

0 

J+ 

0 

J_ 

r_ 

0 

-T+ 

0 

0 

0 

oj — h 

0 

J_ 

0 

J+ 

0 

0 

-T_ 

0 

T+ 

0 

0 

0 

—uj -f h 

0 

-J_ 

0 

-J+ 

0 

0 

0 

0 

-J+ 

0 

J_ 

0 

issi 

-issi + i 

0 

0 

0 

0 

0 

0 

0 

J+ 

0 

-J- 

—issi — i 

issi 

0 

0 

0 

0 

0 

0 

-J_ 

0 

J+ 

0 

0 

0 

issi 

-iSsi + i 

0 

0 

0 

0 

0 

J_ 

0 

-J+ 

0 

0 


ISsl 


we also rescaled frequencies and the magnetic field with 2vf (alternatively, one can think we took wp = 1/2). 

After finding the Green’s functions of interest we replace T± i— Jer/2 and J± i—(Jra ± Al) /2, where in the 
second one (in order to facilitate the comparison with previous results from the literature) we reintroduced the possible 
asymmetry between right and left leads, which stretching our notation is given by AT = (j/^ — jf^)/2. Combining 
these results one gets the following expression for the current: 


^+oo 


I = 


[{J^ + K'^) + dfs) + [Sl (w) - Si (w)] duj 


OJ^ + 


J4 + (J2 + K^ f - 2h^ + [J^ (J^ + K^) + h^f 


(24) 


This expression is completely equivalent to the one S&H reported in their original work. [See Eqs. (7) and 
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(8) of Ref. 20 and match their notation to ours according 
to Fa = Fi + where Fi = J^j., and F(, = In 

particular, it is interesting to consider the case of zero 
magnetic held. One hnds 


^+oo 


I = 


Jl{Jl + K^) 
lo uj^ + {Ji + K^y 


[S; (w) - Si (w)] 


duj 

27T 


(25) 


which is a rather peculiar result, since the dependence on 
Jra has completely dropped out from the problem. This 
integral is elementary in the zero-temperature limit. 


«,= 0 O 


«, = 0 


lim 1 = 

emp —^0 


1 


n, = 1 



r + , 





i r) fvO? 

Jo Cc2 + (J2+A2)2 

sl W 

a 

. 

^ + 

a 


and it can also be carried out at hnite temperature in 
terms of digamma functions [34]. 

Let us focus back on the lead-symmetric case with 
K = 0. When there is no magnetic Held, there is no 
dependence on J^a in the conventional result and it ap¬ 
plies to the case of Jra = 0 in particular. But in that case 
the structure of resembles closely that of the tunnel¬ 
ing term in a simple junction [1], with playing the 
role that was played there by hf and replacing 'tph. 
Therefore, in this limit we know (based on our experience 
with the junction) how to take the h’s into account. And 
we thus know one should expect the consistent result to 
be rather different from the conventional one, except per¬ 
haps for small Jer- Indeed, due to the presence of the n 
factors in the consistent approach, the structure of the 
problem is that of a resonant level attached to a single 
lead. Because of the alternation of different h’s, there 
are no Majorana-like operators that could contribute to 
the transport and give a nonzero current (cf. Ref. 35), 
unlike in the conventional approach. In other words, the 
expected consistent-approach result when Jj-a = 0 is sim¬ 
ply dl/dV = 0 for any value of Jer- It follows that we 
need to treat the h’s consistently for all values of Jra, but 
that requires additional insights. 


FIG. 2. Graphical representation of the states that contribute 
to transport. Intralead processes are indicated by the vertical 
arrows (J±) that conserve h;, while interlead processes are 
given by the horizontal arrows {T±) that conserve hsi- The 
circles indicate the state labels for the single-particle sector of 
an equivalent Gaussian problem that shares the exact same 
processes (see the discussion in the text). 


by processes in and, in particular, do not contribute 
to transport. That leaves only four states as depicted in 
the Hgure. 

From our study of the simple junction [1], we know 
that the physical content of the fi factors is actually to 
avoid contractions between normal and anomalous terms; 
in the present case, however, those are allowed but when 
and only when and Jer processes alternate. By con¬ 
sidering processes at different orders of perturbation in 
one can conclude that we can achieve the same set 
of processes, and also avoid the presence of anomalous 
terms, by modifying the anomalous terms in Eq. (20) 
according to the following prescription: 

dt ^ /t 
ipj —>• -'i’si 
-iph —^ 


B. Consistent Approach 

Let us turn again to the problem of Hnding the nec¬ 
essary expectation values, fully dressed by as in 
Eq. (20), needed to compute the current as in Eq. (23). 
The challenge is to treat the factors and consis¬ 
tently. This is a nontrivial problem, but we are neverthe¬ 
less able to provide an ad hoc solution. In order to achieve 
that, we start by studying the local Hilbert space at the 
impurity site and the structure of the processes that take 
place; this is graphically summarized in Fig. 2. 

Let us consider the possible sets of eigen-expectation- 
values of hd^i^si = 0,1 (cf. Ref. 1). There are eight combi¬ 
nations in total, but those in which they add up to an odd 
number constitute isolated states that are not connected 


while not adjusting those same Helds (only renaming 
them by adding twiddles) in the regular terms, and re¬ 
moving all the fi factors everywhere. In particular, one 
can check processes to fourth order in perturbation, when 
there are combinations in which all possible vertexes en¬ 
ter and the state of the system goes around full circle; 
cf. Fig. 2. We veriHed all these processes are in one-to- 
one correspondence in both formulations. This mapping 
is also argued for in a diHerent way in the Appendix, 
by studying some exactly solvable limits of the two-lead 
Kondo model. 

The mapping has to be applied to a version of 
previous to that in Eq. (20), that is already debosonized 
but in which the factors of still appear explicitly 

in the interlead terms and can be removed ulteriorly. We 
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are thus lead to consider 

(0) + ^L ^^)d) + 

+ (0) + ^/(0)/)- 
-'%{dHi (0)+^t(o)j)_ 
-^(/^^.;(0)+4(0)/) , (27) 


difference with Eq. (20): after applying the prescribed 
mapping, and even though we kept similar notations for 
the fields, there is no connection left to the physical sec¬ 
tors of the theory. 

We also want to consider the local magnetic-field term 
which is naturally rewritten in a symmetric way (see the 
Appendix): 

infield = -h (d^ + pf - l) . (28) 

The current is given by 


where the ± superscripts are used simply to keep track 
of what was the corresponding superscript in the now- 
absent nf ’s of the different terms, but in the calculations 
the two couplings will be taken as having equal numer¬ 
ical values. The factors were removed from the 

couplings by using again a time-dependent gauge trans¬ 
formation. The choice of gauge transformation is this 
time not unique. We adopted the following symmetric 
choice: fii = fj,f = —eE/2 and ^si = = eE/2. Another 

choice could have been fj,i = fj,f = —eV and the other 
two zero; we checked that these and other choices are 
all equivalent. We need to stress that we introduced the 
'ip notation to emphasize an important interpretational 




- (Ml)) + 


(29) 




The resulting calculation is straightforward. Using the 
spinor basis 

= d- d+ f- f+)^ (30) 

and adopting the same notations as for the conventional 
calculation, the local inverse Green’s function for the 
junction is given by 
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where it should be noticed that the voltage now enters explicitly also in the resonant-level-like diagonal-block action 
matrix elements corresponding to d^ and f7 The voltage was written here absorbing a factor of e/2 only to keep the 
matrix expression short; it will be reinserted below. The final expression for the current is 


r+oo 


I = 


cko 


2 (J2 - Jl) {OJ + hf 


27r 


{u; + hf-Jl-{eV/2f +4Jg(uj + h) 


[s; (w) - Ssl (w)] , 


(31) 


where Jss/ds = iPa^Jer)/^ ^''6 the sum and the difference of the squares of the couplings, respectively. Notice that only 
the squares of the couplings enter into the final expression, which means that all the n’s would have appeared squared 
as well. As anticipated, the current vanishes if either J^a or Jer is zero. The integral can be done in general, but it is 
simpler and more illuminating in the zero-temperature limit. Using d (s; — Ssi) /dV = (e/2) d (s; — Ssi) /d (eU/2) —>■ 
e[5 {uj -\- eV 12) -I- 5 (w — eU/2)], we directly write down an expression for the zero-temperature differential conduc¬ 
tance. 


dl 

dU 


= E 


2(e/2^) ( 4 - {h + seVl2f 


S — ± 


+ 


1 [{h + s eV/2f - - {eV/2f] + 4 J2 {h + s eV/2f 

4e (eU/2) (4, - + hf - - {eV/2f] {u: + h) 


27r 


n2 


[Ud + hf -Jl-(eV/2Y +AJl{uj + h) 
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FIG. 3. Comparison of the differential conductance for the 
two-lead Kondo junction, G — dl/dV (in units of the single¬ 
channel conductance quantum, Gq = e^/h), calculated ac¬ 
cording with the conventional procedure or using our consis¬ 
tent scheme. The plot is at zero applied voltage and finite 
magnetic field {h = wp). Notice the convention used for the 
horizontal axis in order to cover the full range of J (see also 
the corresponding plot for the simple junction [1]). Here, dif¬ 
ferently from previous sections, J stands for either Jer or J^a 
while the other (Jra or Jer, respectively) is set to 1; the re¬ 
sulting plots are identical for both cases at zero voltage. The 
dashed line gives for comparison the conventional result di¬ 
vided by a factor of 2. 


C. Comparison of Results 

In what follows we will illustrate the differences be¬ 
tween the results obtained using the conventional and 
consistent approaches discussed in the previous two sec¬ 
tions. We will do so mainly by considering the behavior 
of the differential conductance, first as a function of Jer 
and Jra for specific values of h and eV as shown in Fig. 3 
and then by looking at its behavior as a function of h and 
eV for specific values of Jer and Jra as shown in Figs. 4 
to 6. The physical interpretation and plausibility of the 
consistent-approach results will be presented as well. 

Figure 3 shows the values of the differential conduc¬ 
tance for voltage V = 0 and field h = vp a.s a function 
of Jer while fixing Jra = 1, or vice versa. We denote Jer 
and Jra by J in the two cases, respectively. The two cases 
produce completely identical plots at zero voltage, as can 
be seen from the symmetry of Eqs. (24) and (31). From 
these expressions, one can also see that for the current 
and differential conductance derived using the consistent 
scheme, the symmetry under exchange of Jer and Jra con¬ 
tinues to be present also at finite voltage. On the other 
hand, the symmetry is absent when one puts h = 0 for 
any value of voltage in the conventional expression and 
further the transport is dictated solely by the value of 
Jer. It is to be noted that the order of limits for Jer, Jra, 
V, and h going to zero is in general important and needs 
to be treated carefully, as can be seen by a direct study 


of the expressions for the current. 

From a diagrammatic point of view, one would have 
expected the conventional and consistent plots in Fig. 3 
to be asymptotically equivalent for small J. An expan¬ 
sion in the couplings of the zero-voltage dl/dV derived 
from Eqs. (24) and (31) shows that, unlike the case of 
a junction [1] or the Ising limit considered in the Ap¬ 
pendix, there is no contribution to the current in the 
lowest order. That is the order at which the two calcu¬ 
lations would have matched (as indeed happens for the 
simple junction and the Ising limit). The next order is 
which is the first nonvanishing order for the 
differential conductance and the first one consistent with 
the emerging Jer O Jra duality of the problem at the 
Toulouse point that was discussed above. At this order, 
the conventional calculation is larger than the consistent 
one by a factor of 2; which is the same factor that one 
finds for the (same) next-to-lowest order of the expan¬ 
sion of the consistent calculation as compared with the 
exact direct results in the cases of both the junction and 
the Ising limit. Interestingly, as can be seen from the 
comparison of the dl/dV results for Jer = Jra, the two 
calculational schemes yield, at zero voltage, results that 
are identical except for the aforementioned factor of 2. 
Motivated by these two observations, we also provided in 
Fig. 3 the conventional result scaled by a factor of 1/2. 

Although the tails do not match in Fig. 3, the two 
ways of calculating give the same result as J (i.e., Jer 
or Jra) goes to zero while the other one (i.e., Jra or Jer, 
respectively) is nonzero. However, similarly to the case 
of a simple junction [1], the results are different for fi¬ 
nite J, the difference being most marked for large J. 
While the consistent conductance first increases and then 
falls down, eventually going to zero as J —>■ oo, the con¬ 
ventional conductance continues to grow with increasing 
J and approaches the value 2e^//i. 

One can gain some additional insights by appealing 
to physical arguments. The current is obviously zero if 
Jer = 0, but also if Jer —OO due to the formation of a 
resonating-tunneling state that will block other electrons 
from approaching the junction (the same as what hap¬ 
pens for a simple junction [1], since those arguments are 
not affected by the fact that tunneling now involves also 
spin flip). Therefore, a nonzero current requires the inter¬ 
lead exchange coupling to take some intermediate value. 
The physical picture is similar but more involved con¬ 
cerning the intralead couplings. If Jra = 0, for instance, 
then the parallel Kondo terms which reach their maxi¬ 
mal value at the Toulouse point (measured in terms of 
phase shifts) promote the formation of a strongly bound 
“static doublet” between an electron at each side and 
the impurity. Such a bound state “sits” at the location 
of the junction and Pauli-blocks the passage of a current: 
since electrons tunnel with spin flip, there is always one 
bound-state electron already occupying the site with the 
same spin projection that the tunneling electron would 
have either before or after tunneling. If Jra —>■ oo, in the 
opposite limit, then a complicated resonating-exchange 
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FIG. 4. Contour color maps of the differential conductance, G = dl/dV, in units of (twice) the single-channel conductance 
quantum as a function of applied field and voltage. The two panels compare the conventional calculation (left panel, in units of 
2 Gq) and our consistent calculation (right panel, in units of Gq) for the case of Jer = 1 = Jra- While, conventionally, applying 
a field or a finite voltage have the same effect, that is not the case in a consistent calculation. 


state would form instead, but again the current would be 
blocked because other electrons would be blocked from 
approaching the junction (similarly as for Jgr —> oo). As 
a result, a steady-state current requires intermediate val¬ 
ues of both inter- and intralead couplings. 

In all the limits that involve strong coupling, the in¬ 
tuitive physical picture discussed above does not agree 
with the results of the conventional calculation, but it 
does with the consistent calculation in which the h’s are 
properly taken into account. This is a strong validation 
for the need to manipulate and debosonize models con¬ 
sistently after the initial bosonization. Additionally, it 
should also be remarked that the consistent conductance 
never exceeds the value of one single-channel quantum 
of conductance {Gq = jh). This is also an emergent 
property at the nonequilibrium Toulouse point in addi¬ 
tion to the Jer Jra transport duality, both of which are 
indicated by physical arguments along the lines discussed 
above. 

It is interesting to take a more systematic look at 
the variation of the differential conductance with applied 
field and voltage. In Fig. 4 we show contour color maps of 
dl/dV for the conventional and consistent calculational 
schemes, while fixing the Kondo couplings to the sym¬ 
metric choice Jer = 1 = Jra- For additional clarity, in 
Fig. 5 we show also two half-plane cuts of the differential 
conductance maps: a vertical one at constant magnetic 
field, h/vF = 5, and a horizontal one at constant voltage, 
eV/vp = 5. For these cuts, the horizontal axis is defined 
in terms of either x = eV/h (dashed lines) or a; = h/eV 
(solid lines), respectively. As is clearly evident from the 
figures, the effects of field and voltage turn out to be ex¬ 


actly equivalent in the conventional calculation (as seen 
from the 90° rotational symmetry of the first contour plot 
or from the complete overlap of the two “conventional” 
traces in the cuts). However, the consistent calculation 
yields inequivalent dependencies on field and voltage. 

Both ways of calculating show a splitting of the zero- 
bias differential-conductance anomaly due to the finite 
magnetic field when plotting as a function of voltage, but 
the peaks are sharper and more asymmetric in the consis¬ 
tent calculation. The contrast between the two results is 
even greater at finite applied voltage when plotting as a 
function of magnetic field. In the consistent case and for 
low applied voltages, only part of the zero-bias anomaly 
splits, while a relatively broad relic of it remains pinned 
at zero bias. As a result the contour plots are star- or 
butterfly-shaped, instead of being cross-shaped as in the 
conventional calculation. For small x (as compared with 
Jer and Jra)j the differential conductance can addition¬ 
ally be described as showing a small ‘deep’ developing 
right after/before the ‘peak’ when plotting as a function 
of voltage or field, respectively. In other words, besides 
the ‘ridges’ in the contour maps for \h\ = 0, \eV\, there 
are also ‘furrows’ at |2h| = \eV\ which are absent from 
the conventional calculation. This can be seen analyti¬ 
cally by looking at the small-coupling expansion of the 
differential conductance. 

To take a closer look at the central part of the contour 
maps, we plot in Fig. 6 horizontal and vertical cuts going 
through h = V = 0 (solid and dashed lines, respectively). 
In the conventional calculation, the differential conduc¬ 
tance reaches the maximum value of 2Gq at zero field 
and bias, and decreases towards zero as either of them 
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FIG. 5. Differential conductance for the two-lead Kondo junc¬ 
tion, G = dl/dV (in units of the single-channel conductance 
quantum), calculated according with the conventional proce¬ 
dure (in red) or using our consistent scheme (in blue). The 
plots are either at a hnite constant field {h = 5uf) and as 
a function of applied voltage (dashed lines), or at a finite 
constant applied voltage [eV = 5wf) and as a function of 
magnetic field (solid lines). For the conventional calculation, 
both plots are identical and the dashed line is covered by the 
solid one. Notice the convention used for the horizontal axis 
in order to cover the full range of x (defined in terms of either 
X = eV/h or X = hleV, with the field that is kept constant 
being the one in the denominator). 


increases (once again, the symmetry of the behavior with 
h and V make the solid and dashed lines coincide with 
each other). On the other hand, the consistently calcu¬ 
lated differential conductance reaches a maximum value 
of Gq at the origin—only half as tall—and it decreases 
towards zero differently with field or with applied voltage, 
with the latter being the slowest decay and the only one 
of these four curves not following a steepest descent. No¬ 
tice in addition, from the contour maps, that only along 
the diagonal ridges would the differential conductance 
not go to zero asymptotically. Let us point out that the 
way the asymptotic values are approached in all cases is 
as a power law instead of the expected logarithmic tail 
[23, 36]. This is a peculiarity of the Toulouse limit, in 
which the voltage or local magnetic field can never be 
larger than the parallel Kondo couplings and the band¬ 
width is infinite. So the standard argument in which 
such energy scales stop the renormalization-group flow of 
the couplings does not apply. Moreover, already in the 
conventional framework and for equilibrium situations, a 
standardly formulated renormalization-group scheme [37] 
is not compatible with the nature of the Toulouse-limit 
fixed-point manifold, which calls for a careful reformula¬ 
tion [29, 38]. 

Let us also stress and comment on the differences be¬ 
tween the ways in which applied voltage and temperature 
enter in both calculations. In the conventional calcula¬ 
tion, both enter into the expression for the charge cur¬ 



eF [2vp] or h [2vp] 

FIG. 6. Differential conductance for the two-lead Kondo junc¬ 
tion, G = dl/dV (in units of the single-channel conductance 
quantum), calculated according with the conventional proce¬ 
dure (in red) or using our consistent scheme (in blue). The 
plots are either at a zero field and as a function of applied 
voltage (dashed lines), or at a zero voltage and as a function 
of magnetic field (solid lines). For the conventional calcula¬ 
tion, both plots are identical and the dashed line is covered 
by the solid one. 


rent, Eq. (24), through the thermal factors si (w) and 
Ssi (w) only, but not via the kernel multiplying them 
(the spectral function of a Majorana fermion in this 
case). The latter is independent of both eV and Temp, 
which is ascribed to the quadratic nature of the prob¬ 
lem at the solvable point [34]. Only as one moves away 
from the solvable point, the voltage and temperature 
will explicitly enter into the kernel of the integrand [39]. 
This is in contradistinction to the consistent calcula¬ 
tion in which the voltage enters explicitly into the ker¬ 
nel of Eq. (31), but not the temperature. Thus, not 
only voltage and magnetic field, but also voltage and 
temperature, interplay differently as compared with the 
previously accepted, conventional results for the non¬ 
equilibrium Toulouse point. If one were to imagine that 
these kernels capture the Kondo-resonance part of the 
spectral function of some parent Anderson-type impu¬ 
rity model, one could see how the resonance reacts to 
applied field or voltage. On the one hand, in both cal¬ 
culations the resonance would shift from zero frequency 
with magnetic field. On the other hand, the resonance 
would not depend on voltage in the conventional picture, 
but it would split with the applied drain-source bias in 
the consistent picture. The latter would be in agreement 
with experimental results and the accepted phenomeno¬ 
logical picture of Kondo transport out of equilibrium [40] 
(and as first predicted in Ref. 41). 
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V. CONCLUSION AND FINAL REMARKS 

In the companion work of Ref. 1, we presented a consis¬ 
tent bosonization-debosonization program in which we in¬ 
troduced the n factors, defined here in Eq. (9), to assist in 
making the results consistent after performing transfor¬ 
mations in the bosonic language and later debosonizing 
models that include terms with single-point non-normal- 
ordered operators. How these factors should be treated 
depends on the physical setup being considered and needs 
to be studied on a case by case basis. Sometimes the con¬ 
ventional way (i.e., h —)■ 1 for all of them) could be the 
one consistent with the problem, but other times the con¬ 
sistent treatment is different. Moreover, this treatment 
can be connected to the choice of boundary conditions 
that is dictated by the problem (cf. Ref. 1). 

We applied these ideas here to the important case of 
quantum impurity problems. In particular, we focused 
our attention on the two-lead Kondo model of a junction 
out of equilibrium. By considering certain regimes of the 
problem, we were able to argue that the conventional 
way of calculating does not produce consistent results, 
while a different treatment of the n factors seems to fix 
those problems (as it did for the case of the simple junc¬ 
tion problem [1]). Moreover, the calculations can then be 
carried out in the full regime of parameters of the system. 
This way, the key insights of the work by S&H (and also 
by E&K) can be retained and the calculations hxed to 
produce consistent results (our method of solution was 
ad hoc and there is no exact solution to refer to as in 
the case of the junction, but we do know it interpolates 
between consistent limits). We thus were able to make 
a number of predictions for the transport characteristics 
of the two-lead Kondo model that can, in principle, be 
looked for in experiments. 

Certainly, more work should be done along this line 
following the developments of the literature of the past 
two decades, and we are already exploring some direc¬ 
tions. To name but a few: (i) the consistent solution 
can be explored further, including additional aspects of 
transport (such as the noise spectrum and possibly ther¬ 
mal transport), alternatives such as charge sensing [42], 
and also the thermodynamics; (ii) one can study ac- 
drive effects that are realizable in experiments [43, 44]; 
(iii) a study of multiterminal models would provide new 
insights, cf. Refs. 31, 45; (iv) perturbation around the 
solvable Toulouse point needs to be considered anew, 
cf. Refs. 39, 46; (v) connections to other approaches such 
as boundary GET can bring in synergy [47, 48], show¬ 
ing for example how to possibly extend those methods 
to nonequilibrium transport problems; (vi) ditto for ap¬ 
proaches that exploit the connections to integrability [49] 
or to renormalization ideas [50]; (vii) one could combine 
our approach with a hnite-size bosonization analysis to 
further bridge with GET and numerical renormalization 
ideas [28]. The list can go on (even though we restricted 
it to quantum impurity problems only). Given the con¬ 
tinued challenge posed by the need to better understand 


strongly correlated quantum systems out of equilibrium, 
the consistent Toulouse-point solution will play an im¬ 
portant role as a reference case for a class of problems in 
which a lot is still not well understood and there is a lack 
of exact results to guide the theoretical developments. 
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Appendix: Some Exactly Solvable Limits 

In the following, we consider two different limits that 
can be rigorously treated via exact direct calculations. 
They further motivate the prescription mapping dis¬ 
cussed in the main text to solve for the transport within 
the consistent approach. 


1. f-axis Ising Limit 

We call the Ising limit of the anisotropic Kondo model 
that in which all the spin-spin exchange interactions take 
place along a single axis. We will choose it to be the (as 
is often done for the transverse-field Ising model). The 
virtue of this limit is that it is exactly solvable: the elec¬ 
trons still interact with the impurity but the latter does 
not have dynamics (i.e., no spin flips take place along the 
X axis). As a result, all one has to do is to solve for each 
possible impurity orientation and average over the re¬ 
sults (this is reminiscent of the treatment of the boundary 
sine-Gordon model that one obtains at the solvable point 
of the problem of a classical impurity in a Tomonaga- 
Luttinger liquid [31, 52]; cf. Ref. 53). When the impurity 
is frozen, the only remaining degrees of freedom are the 
electronic ones and the problem becomes Gaussian and 
thus exactly solvable in a direct way. We shall keep the 5 
axis as the quantization axis for the electrons, as this will 
show some structure that will help us understand how to 
deal with the full two-lead Kondo model in the language 
of Abelian bosonization. 

In the zero-field case, and setting in Eq. (1) all the 
non-a; couplings to zero, = 0, we are left with 
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(the time dependence is implicit) 

H = '^( fnUx + H^) , (A.la) 

= 'ipl^ (x) {-ivFdx) ipcre (x) , (A.lb) 

= Ju' ^imp (0) (0)) , (A.lc) 

where we have made use of the relation 2S';^p5'g[gj. = 

5'fmp {Stiec + 5'elec)- Hermiticity requires = 

Jgr- It is a simple exercise to solve the model directly and 
find a closed expression for the differential conductance 
(for brevity, we do not quote that result here). 

We are interested in bosonizing, changing boson basis, 
and debosonizing in the same way as we did for the full 
model (gauging the voltage in and out from the couplings 
also in the same way). Of course, no unitary transforma¬ 
tion is required this time since S'j^p commutes with the 
Hamiltonian and thus plays a simple spectator role. The 
Kondo part of the Hamiltonian is hnally rewritten as (all 
fields are at time t and x = 0) 

+ i’si’sl) - 

- ^er^^imp (V'lV's + ^l^l) • (A.2) 

One could now set all h —>■ 1 and proceed convention¬ 
ally to calculate the differential conductance. The result 
is that the expression differs from the direct calculation 
in a way that parallels what we discussed for a simple 
junction [1]. In particular, a small-coupling expansion 
shows that the results match to lowest order, C>[(J^,.)^], 
but the conventional result is twice bigger than the di¬ 
rect one at next-leading order, as was the case for the 
simple junction. Therefore, one needs to treat the h’s 
more carefully. 

The practical problem that arises is that the factors 
and nfi introduce complicated dynamics into the Hamil¬ 
tonian because there are linear terms in and 
present as well. From our study of the simple junction [1], 
we know that the physical content of these factors is ac¬ 
tually to avoid contractions between regular and anoma¬ 
lous terms of the same type (i.e., inter- or intralead). 
Here the situation is more complicated because regular 
and anomalous terms of a given type can be combined 
provided there is an intervening term of the other type, 
and vice versa. One can check this conclusion order by 
order via a matching of perturbative expansions as we 
did for the case of the simple junction [1]. 

But since the present problem is directly solvable in 
terms of the original fermions, we know there exists a way 


of organizing the perturbation theory as if the model was 
purely Gaussian. We shall thus focus on the structure of 
the terms while comparing new and original fermions. 
We do that by looking at the four Klein factor rela¬ 
tions [see Eqs. (lOa)-(lOh)] that involve F^: (i) for the 
regular terms, FsiF^ = F^^F^f and FiFj = F^^F^f] 
and (ii) for the anomalous terms, F^^F^ = F^^F^^, and 

Fi F^ = —F^^F^f- The other four are just the Hermitian 
conjugates of these. 

A first observation is that in terms of the original 
fermions there are no anomalous terms present in the 
model; see Eq. (A.l). This prompts us to attempt a 
mapping in which we do not change the regular terms 
but ‘regularize’ the anomalous ones. Notice that while 
F^ appears in all four terms, no single original Klein 
factor appears four times; while ^1?^ repeats in the reg¬ 
ular terms, does in the anomalous ones. A second 
observation is that the pair of spin-down original Klein 
factors repeat in regular and anomalous terms, but they 
exchange roles as to which one goes in the intralead pro¬ 
cess and which in the interlead process in each case. 

These two observations indicate that we can achieve 
the same perturbative processes (and, as a bonus, avoid 
the presence of anomalous terms) by modifying only the 
anomalous terms according to the following prescription: 

^ 

'tpj —-'ipsi (A.3) 

iph 

and removing all the n factors. To keep the same phase 
conventions as we used for the matching of Klein factor 
bilinears, we also need to introduce a minus sign in all the 
terms (this is not essential as the current is not sensitive 
to it). Explicitly, one has 

Hk = 2'^R‘^iW + 

+ 2'^L'S'knp + 

+ + 

+ + ^l^i) (A.4) 

and, trivially since the ^’s can be put in one-to-one cor¬ 
respondence with the original fermions in Eq. (A.l), all 
the results will be the same as in the direct solution. We 
stress that, even though we kept the notation with s, I, 
and s/, these are different fermions and there is no direct 
connection left to the physical sectors of the theory (we 
introduced the Tp> notation to emphasize this point). 

2. Flat-band Limit 

If one introduces a lattice discretization, the flat-band 
limit is the limit of zero hopping between sites. Then all 
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sites in each band are independent fermionic degrees of 
freedom, except for the sites at a;o, the location of the 
impurity, which are connected by the tunneling terms 
in the Hamiltonian. We could consider this limit for the 
Toulouse-point Hamiltonian with its Kondo part given by 
Eq. (20) or for the BdB-mapped cc-axis-Ising Hamiltonian 
corresponding to Eq. (A.2). Both yield equivalent models 
in the flat-band limit; for brevity, we frame our presen¬ 
tation around the second case. Taking the band energies 
to be zero, we introduce the notation ipl (xq) —> cj, for 
i' = s,l, si (or d'l’ —>■ cl in the first case) to emphasize 
the lattice nature of the problem and the sector of the 
Hamiltonian containing x = xq turns into the follow¬ 
ing three-site model (taking the expectation value of the 
“spectator” impurity spin and absorbing it in redefined 
coupling constants), 

Hss = -k -k -k 

+ (44 + QCs) - Jl^n~ (cjcs -k . 

(A.5) 


Since now we are dealing with lattice fermions, we can 
directly identify n^. 

This model has an eight-state Hilbert space split 
into two particle-number-parity sectors. The states 
with even number of particles are all disconnected 
degenerate states with zero energy. The states 
with odd particle number hybridize and can be 
diagonalized (using exact diagonalization) into four 


eigenstates with energies ±^—b/2 ± ^Jlp- jA, — c where 


h = - 


m m + (4"r) + 


and 


r /a rr _ 7a 7r 
V^ra'^ra '^er'^er/ • 

Applying the prescribed changes for the anomalous 
terms, this Hamiltonian turns into the equivalent of 
Eq. (A.4), which is now a four-site model with no anoma¬ 
lous terms. 


= - JA 
-Jt 


[clh + - Jra (4cs/ + 4/^^) “ 

(c%i + chc^ - Jl, (cjcs -k 4 q) ■ 


(A.6) 


This model conserves particle number and, in the single¬ 
particle sector, its spectrum coincides with the odd sec¬ 
tor of iJss (cf. Fig. 2). Moreover, the same conclusion 
remains true when we introduce a local magnetic field as 
in Eqs. (19) or (28), appropriately rewritten. This serves 
as a check of how to correctly normalize the local-field 
term after the mapping. In other words, the Hamiltonian 
matrices in those two sectors are unitary equivalent, 

^’odd^^Ss^’odd = PlpH/^gPlp ■ (A.7) 

Due to the absence of interactions in i74s, its many- 
particle physics can be calculated in terms of a Green’s 
function for a single fermion injected into an empty band 
[54], and thus the single-particle sector is the crucial one 
to determine the full dynamics. Notice, in addition, that 
the way this kind of Green’s functions at the bare level 
enter into the transport calculations of the main text is 
via their inverses. As such, the propagator of any hy¬ 
bridized local or flat-band degree of freedom will not 
require independent regularization and will inherit its 
causal properties (i.e., its Keldysh structure [33]) from 
other extended-band degrees of freedom; and thus the 
difference between “empty-band” and “degenerate-gas” 
Green’s functions does not enter the calculations. Thus, 
for those physical properties whose calculation requires 
only the two-point Green’s functions and, by extension, 
only the single-particle sector, the Hamiltonian mapped 
via the prescription should work fine. And in the partic¬ 
ular example above we were able to show explicitly that 
the single-particle sector of the mapped model reproduces 
the nontrivial part of the original spectrum. 

In summary, the mapping we prescribed for the anoma¬ 
lous terms works in both cases, namely the linear- and 
the flat-band limits. In the case of the two-lead-Kondo- 
model Toulouse point, we encountered two linear-band 
and one local degrees of freedom, which is a combination 
of these two cases we just discussed. One can resort to 
neither bosonization nor exact diagonalization to prove 
the mapping rigorously, but it is nevertheless justified on 
physical grounds as a combination of these two limits and 
as argued also in terms of processes in the main text. 
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